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Abstract 



We consider the U{n) x U{m) symmetric Lagrangian to describe the 
finite-temperature phase transition in QCD in the hmit of vanishing quark 
masses with n = m = Nj flavors and unbroken anomaly at T^. We compute 
the Renormalization Group functions to five-loop order in Minimal Subtrac- 
tion scheme. Such higher order functions allow to describe accurately the 
three-dimensional fixed-point structure in the plane (n, m), and to reconstruct 
the line n^(m, d) which limits the region of second-order phase transitions by 
an expansion in e = 4 — d. We always find n"'"(m, 3) > m, thus no three- 
dimensional stable fixed point exists for n = m and the finite temperature 
transition in light QCD should be first-order. This result is confirmed by the 
pseudo-e analysis of massive six-loop three dimensional series. 

PACS Numbers: 12.38.Aw; ll.lO.Kk; 64.60.Fr 
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Typeset using REVT^ 



I. INTRODUCTION 



The phase diagram of QCD is characterized by a low temperature hadronic phase with 
broken chiral symmetry and an high temperature phase with deconfined quarks and gluons, 
in which chiral symmetry is restored. The nature of the transition between these two phases 
depends on the QCD parameters, as the number of flavors and quark masses. In the limit 
of zero quark masses such phase transition is essentially related to the restoring of chiral 
symmetry (see e.g. the reviews [1]). 

The QCD Lagrangian with Nf masslcss quarks is classically invariant under the global 
flavor symmetry U{1)a x SU{Nf ) x SU{Nf ) [2]. The axial U{1)a symmetry may be broken 
by the anomaly at the quantum level, reducing the relevant symmetry to SU (Nf) x SU (Nf) x 
Z{Nf)A [2]. At T = the symmetry is spontaneously broken to SU{Nf)v with a nonzero 
quark condensate. With increasing T, a phase transition characterized by the restoring of 
the chiral symmetry is expected at Tc- To parameterize this phase transition a complex 
Nf-hj-Nf matrix is introduced as an order parameter. The most general renormalizable 
three-dimensional U{Nf) x U{Nf) symmetric Lagrangian is [2,3] 

CuiNf) = Tr(a^<l>t)(9^$) + rTr$t$ + ^ (Tr$t$)' + ^Tr ($ t$)' , (1) 

which describes the QCD symmetry breaking pattern only if Vq > [4]. 
If the anomaly is broken at T^., the effective Lagrangian is [2,3] 

J^su{Nf) = j0.u(Nf) + wo(det $ + det $^). (2) 

The effect of non-vanishing quark masses can be accounted for by adding a linear term 
rriij^ij in the Lagrangian [1-3], that acts as a magnetic fleld in a spin model. 

The mean-field analysis of the U{Nf) x U{Nf) Lagrangian Eq. (1) predicts a second- 
order phase transition everywhere the stability conditions vq > and NfUo + vq > are 
satisfied. However, according to Renormalization Group (RG) theory the critical behavior 
at a continuous phase transition is described by the stable fixed point (FP) of the theory [5]. 
The absence of a stable FP indicates that the transition cannot be continuous, even though 
mean-field suggests it. Therefore the transition is expected to be first-order (see e.g. [6]). 

The U{Nf) x U{Nf) Lagrangian was studied at one-loop in e — A — d expansion [2,7,8] 
and at six-loop in the massive zero-momentum renormalization scheme directly in ci = 3 [4] . 
No stable FP was found for all values of Nj > 2, concluding for a first-order phase transition. 
Anyway both the used approximations have intrinsic limits. The one-loop e expansion, as 
discussed in Ref. [4], provides useful qualitative indications for the description of the RG 
fiow, but it fails in describing quantitatively the right three-dimensional behavior. The 
major drawback of fixed dimension expansion is that the numerical resummation techniques 
necessary to extract quantitative informations allow to explore a large, but limited, region 
in the space of coupling constants (e.g. in Ref. [4] the resummation results to be effective 
only in the region —2 < u,v < 4, where u,v are the couplings used in Ref [4]). One 
cannot exclude a priori that a FP may be outside the accessible region of effectiveness of 
resummation in = 3. These problems are absent in e expansion since no resummation is 
needed to find the FP's, being series in e. 
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For these reasons, we extend the e expansion series to five loops. We consider the 
U{n) X U{m) generalization of Lagrangian (1), where $ is a n-by-m complex matrix. To 
understand why we decide to study this more general model let us consider the already 
known one- loop (5 functions [7]: 

v) = — + (nm + 4)m^ + 2(n + m)MV + 3f ^ , P^iu^v) — —ev -\- Q>uv -\- {n -\- w^jv"^ . (3) 

For n = m a couple of FP's with non- vanishing and negative v exists only for n < \/3 + 0(e) 
[2], suggesting a first-order phase transition for those systems, as light QCD, having vo > 0. 
Anyway, if one considers the model with n ^ m a. more complicated structure of FP's 
emerges. Two FP's, the Gaussian one {u* = v* = 0) and the 0{2nm) one {v* = 0), always 
exist. For n > n~^{m, d) and n < n^{m, d) other two FP's appear which we call [/"*" and U~ , 
whose coordinates at one-loop read 

"'± — in ^1 "± — D ' \ ) 

with 

Bmn — nm^ — hn-\- mn^ — 5m , Amn = 36 — — 2mn — r? , (5) 
R^n = 24: + m^ - lOmn + , Dmn = 108 - 8m^ - 16mn + m^n - Sri^ + 2m^n^ + mn^ . 

The V* coordinates of t/^ for n > are positive, as it should be to provide the right 
symmetry breaking of QCD. Requiring Rmn > 0, we have = 5m± 2\/6\/m^ — l + 0{e)^. 

The stability properties of these FP's at fixed m (for the physically relevant case m > 1) 
are characterized by the following four different regimes (note the analogy with the 0{n) x 
0{m) model [9]): 

1) For n > n^{m,d), there are four FP's, and C/+ is the only stable. Both have 
V* > 0. 

2) For n~{m,d) < n < n'^{m,d), only the Gaussian and the Heiscnbcrg 0(2nm)- 
symmetric FP's are present, and none of them is stable. Thus the transition is expected 
to be first-order. 

3) For nii-(m, (i) < n < n~{m,d,), there are again four FP's, and is the stable one. 
However at small e for m < 5, it has v* < 0. For vq > a first-order transition is 
expected. For this reason we will never consider the value of n^{m, d). 

4) For n < nnim, d), ^ there are again four FP's, and the Heisenberg 0(2mn) -symmetric 
one is stable. 



^Note that it is possible to study the generalized U{n) x U{m) model at fixed m in a 1/n expansion, 
since it has a FP for n = oo, contrarily to the U{n) x U{n). 

^The value of nu{m) may be inferred from Refs. [10,11], where it was shown, on the basis that at 
a global 0{N) FP all the spin four operators have the same scaling dimension, that the 0{N) FP 
is stable only if iV < AT^, with A^^ ~ 2.9 [6,11,12]. Thus the 0{2nm) FP is stable for n < 1.45/m 
and it is not expected to play any role at the QCD phase transition. 
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Now it is possible (and one has to check!) that the actual value of n+(m, 3) is lower 
than m, providing a stable FP and consequently a new universality class for U {Nf) x U (Nf) 
symmetric models. To give a definitive answer to this question, high order calculations are 
required, since low order ones lead to erroneous conclusions, as we shall see. However we 
anticipate that we do not find any stable FP, supporting the results of Rcf. [4]. 

The paper is organized as follows. The U{n) x U{m) model is analyzed at five-loop in 
e expansion in Sec. II. In Sec. Ill, the model is reanalyzed with pseudo-e expansion at 
six-loop order in massive zero-momentum renormalization scheme. Sec. IV summarizes our 
main results. In the appendix we briefly discuss the effect of the anomaly. 



II. FIVE-LOOP e-EXPANSION OF U{N) x U{M) MODEL. 

We extend the one-loop e expansion of Refs. [7] for the RG functions of the U (n) x U (m) 
symmetric theory to five-loop. For this purpose, we consider the minimal subtraction (MS) 
renormalization scheme [5] for the massless theory. We compute the divergent part of the 
irreducible two-point functions of the field of the two-point correlation functions with 
insertions of the quadratic operators <l>^, and of the two independent four-point correlation 
functions. The diagrams contributing to this calculation are 162 for the four-point functions 
and 26 for the two-point one. We handle them with a symbolic manipulation program, which 
generates the diagrams and computes the symmetry and group factors of each of them. We 
use the results of Ref. [13], where the primitive divergent parts of all integrals appearing in 
our computation are reported. We determine the renormalization constant Z$ associated 
with the fields $, the renormalization constant Zf of the quadratic operator <1>^, and the 
renormalized quartic couplings u,v. The functions ?7<^ and rjt are determined using 

the relations 



uo,vo 



(6) 



uo,vo 



d log Z$ 



dlogji 



Uo,Vo 



Vtiu,v) 



dlogZt 



dlog n 



(7) 



Uo,Vo 



The zeroes {u*,v*) of the /3 functions provide the FP's of the theory. In the framework 
of the e expansion, they are obtained as perturbativc expansions in e and then are inserted 
in the RG functions to determine the e expansion of the critical exponents. 



rj ^r)^{u*,v*), V = {2- r]^{u\v*) - r]t{u\ v*)) \ 



(8) 



A. RG functions 



The five-loop expansions of the {3 functions are given by 
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Pu — —^u + {nm + A)u^ + 2(n + m)uv + 3f ^ — |(7 + 3mn)u'^ — ll(m + n)u^v 

_41^^^2 _ ^ ^^^3 ^ ( 740+461mn+33m^n^ ^ C(3)(33 + 15mn)) 

^ j- 659n+79m^n+659m+7W ^ 48((3)(m + Tl)) u\ 

+ ( 2619+1210mn+230n^+230m^+3n^m^ ^ ^8^(3) (7 ^ ^2^2 

+ (x(20n + m^n + 20m + n^m) + 36C(3)(m + n)) wt;^ 

^ j- 425+20m^+153mn+20n^ ^ g^^3^^4 ^ ^4 ^ ^ ^« (g) 

v) = -ev + Quv + {n + m)v'^ - ^^^^u\ - ll(m + n)uv'^ - 3^^v^ 

jj- 821+184mn-13m^n^ ^ 12C(3) (7 + m/l)) ii^ 
^ j- 1591n-35m^n+1591m-35mn^ ^ 72C(3)(m + u)) V^V 
^ ( 211+9m^+73mn+9n^ ^ 24C(3)(4 + mn)) ^Xv' 

^ j^ 295n+13m^n+295m+13mn^ ^ 9C(3)(m + n)) v"^] + (31 + (31 (10) 

The coefficients (3"^ , (3q , , are very long and not really illuminating. We do not report 
them here, but they are available on request to the authors. The same is true for the RG 
functions rj^j, and rjt to five-loop, that we calculated but never used, since we did not find 
evidence for a FP in the space of parameters of interest. 

We have checked that for f = the scries reduce to the existing O(e^) ones for the 
0(2nm)-symmetric theory [14]. For m = 1 and any n (or viceversa, given the symmetry 
under the exchange n ^ m) the theory is equivalent to an 0{2n) in the coupling u + v, so 
the series satisfy the relation Pu{z + y, z — y;n,m = 1)+ (3y{z + y, z — y;n,m = 1) = Po{2n){z), 
where Po{2n){z) is the /3-function of the 0(2n) model [14]. 

B. Estimates of n+(m, 3) 

From the above reported series, the e expansion of n^{m) may be calculated to 0{e^). 
n^(m) may be expanded as 

n±(m) = n^(m) + nf(m) e + n^{m) + nf(m) + n^(m) + O(e^), (11) 

and the coefficients nf{m) are obtained by requiring 

d{(3u, (3v 



Pu{U'*:V*;n^) — 0, P-i,{u*,v*;n^) — 0, and det 



d{u, v) 



{u*,v*;n^)^0. (12) 



For generic values of m, the expression of n^(m,4 — e) is too cumbersome in order to be 
reported here. We only report the numerical expansion of n"*" at fixed m = 2, 3, 4, i.e. 



n 



+ 



n 



+ 



(2, 4 - e) = 18.4853 - 19.8995e + 2.92606^ + 4.6195e^ - 0.7182e^ + O(e^) , 
(3, 4 - e) = 28.8564 - 30.0833e + 6.5566e^ + 3.4056e^ - 0.7958e^ + O(e^) , 



n+(4, 4 - e) = 38.9737 - 40.2386e + 9.6089e^ + 3.0505e'' - 0.6156e* + O(e^) , (13) 
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TABLE I. Estimates of n+ for several m with varying the number of loops. 



m 




l/n+ 






a 




final 




3-l<)C)p 


l-looj") 


5-l()op 


3- loop 


-l-loo]") 


5-loop 




2 


6.01 


4.95 


4.56 


4.98 


4.55 


4.44 


4.5(5) 


3 


9.94 


8.38 


7.76 


8.07 


7.53 


7.40 


7.6(8) 


4 


13.7 


11.6 


10.7 


11.0 


10.3 


10.2 


10.5(1.1) 



In order to give an estimate of n~^{m,3) such series should be evaluated at e=l. A linear 
extrapolation of the two-loop contribution leads to the wrong conclusion that n'^{m, 3) < m, 
i.e. that the transition is continuous. This is the anticipated statement that high-loop 
computation are needed to have a conclusive result. Anyway, a direct sum of the five-loop 
series is not effective, since they are expected to be divergent. The high irregular behavior 
with the number of the loops makes also a Borel-like resummation non effective as well (in 
fact Pade-Borel resummation leads to unstable results). Thus we try to extract from Eqs. 
(13) better behaved series by means of algebraic manipulations. 
This may be done considering (as in Ref. [15]) 

l/n+(2, 4 - e) = 0.0541 + 0.0582e + 0.05416^ + 0.03556^ + 0.0172e^ + O(e^) , 
l/n+(3, 4 - e) = 0.0347 + 0.0361e + 0.02986^ + O.OlSSe^ + 0.0095e^ + O(e^) , 
l/n+(4, 4 - e) = 0.0257 + 0.0265e + 0.0210e^ + 0.0132e^ + 0.0067e^ + C>(e^) , (14) 

whose coefficients decrease rapidly. Setting e = 1 we obtain the results reported in Table I. 

Another method, firstly employed for 0{m) x 0{n) models in Ref. [16], use the knowledge 
of n"'"(m, 2) to constrain the analysis at e = 2, under the (strong) assumption that n"'"(m, d) 
is sufficiently smooth in d at fixed m. n'^{m, 2) may be conjectured further assuming that the 
two-dimensional LGW stable FP is equivalent to that of the NLcr model for all n > 1 except 
n — 1 [5]. Since the NLu model is asymptotically free, wc conclude that n^[m, 2) = 1. The 
knowledge of n~^{m, 2) may be exploited in order to obtain some informations on n'^{m, 3), 
rewriting the perturbative series for n"*" (m, 4 — e) in the following form 

n+(m,4 - e) = 1 (2 - e) a(m,e), (15) 

where 

a(2, e) = 8.743 - 5.579e - 1.3266^ + 1.646e=^ + 0.464e^ + O(e^) , 

a(3, e) = 13.928 - 8.078e - 0.760e^ + 1.323e^ + 0.263e^ + O(e^) , 

a(4, e) = 18.987 - 10.626e - 0.5086^ + 1.271e^ + 0.328e^ + O(e^) , (16) 

whose terms are better behaved than the original series, but not decreasing. We consider the 
series a(2,e)~^ obtaining a more "convergent" expression, which may be estimated simply 
by setting e — 1. The results of this constrained analysis are reported in Table I. Note that, 
although the several not completely justified assumptions we made, the final obtained series 
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TABLE II. Pade table for n"'"(2,3) in pseudo-e expansion. The location of the positive real 
pole closest to the origin is reported in brackets. 





iV = (J 


iV = 1 


A" = 2 


iV = 3 


N = [ N = 5 


M = 





18.485 


5.219 


4.769 


4.596 


4.542 4.527 


M = 


1 


10.762 


4.753[29.5] 


4.487[2.59] 


4.518(3.23] 


4.522(3.72] 


M = 


2 


8.190 


4.579[7.43] 


4.518[3.2i] 


4.523 




M = 


3 


6.921 


4.533[4.99] 


4.522[3.65] 






M = 


4 


6.181 


4.524[4,3oi 








M = 


5 


5.708 











are highly stable with changing the number of the loops. Obviously this is not an evidence 
favoring the validity of the assumptions, but it is a very convincing argument to ensure the 
goodness of our estimates. 

As final estimates we quote an average of the five-loop results (that are quite close) and 
as error bar the maximum difference with the fourth order ones. For all the considered value 
of m we have n'^{m, 3) > m, thus the transition for U{n) x U{n) models is expected to be 
first-order. We also check that n^{m, 3) > m for higher values of m. 

Since the new couple of FP's does not exist for finite temperature transition of light 
QCD, we do not report their expansion in terms of e and the exponents characterizing the 
critical behavior for n > n'^. Anyway they may be obtained from the series we reported and 
from those that are available on request. 

III. PSEUDO-e EXPANSION 

In this section we analyze the six-loop zero-momentum massive three-dimensional series 
with the so-called pseudo-e expansion method [17], since it provided the best results for the 
marginal spin dimensionality in spin models (see e.g. Refs. [15,6] and references therein). 
The P functions for generic n and m were calculated at six-loop in Ref. [4] (but they were 
not reported there). 

The idea behind the pseudo-e expansion is very simple [17]: one has only to multiply 
the linear terms of the two f3 functions by a parameter r, find the the common zeros of the 
/?'s as series in r and analyze the results as in the e expansion. The critical exponents are 
obtained as series in r inserting the FP's expansions in the appropriate RG functions. Note 
that, differently from the e expansion, only the value at r = 1 makes sense. The reason for 
which it works well is twofold: first in the three dimensional approach at least one order 
more in the loop expansion is known, second, and more important, the RG functions are 
better behaved in the massive approach [5,17]. 

Following the recipe explained in the previous section for the MS scheme, we obtain 

n+(2, 3) = 18.4853 - 13.2663r - 0AA99t^ - 0.1735r3 - 0.0537r^ - 0.0144r^ , 
n+(3, 3) = 28.8564 - 20.0555t - 0.3092r2 - 0.2609r^ - 0.1444t^ - 0.0968r^ , 
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n+(4, 3) = 38.9737 - 26.8257r - 0.2690r2 - 0.3582r^ - 0.2199t^ - 0.1526r^ . 



(17) 



At least up to the known order, such expansions do not behave as asymptotic with factorial 
growth of coefficients and alternating signs. So one may apply a simple Padc resummation 
[15]. The results of the [N/M] Fade approximants are displayed in Tab. 11 for m = 2 . 
Several approximants have poles on the positive real axis. Anyway all these poles are "far" 
from r = 1, where the series must be evaluated. Thus one may expect the presence of a pole 
not to influence the result at r = 1. Anyway for the cases m — 3,4: some Fade have poles 
at r < 2, that must be discarded in the average procedure. We choose as final estimate the 
average the six-loop order Fade without poles at r < 2 (excluding those with = 0, giving 
unreliable results), and as error bar we take the maximum deviation from the average of 
four- and five- loop Fade (as in Ref. [15]). Within this procedure we have n+(2, 3) = 4.52(7), 
n+(3,3) = 7.98(25), n+(4,3) = 11.1(4). 

The final results are in very good agreement with those of the previous section from 
a completely different approach. This is a clear evidence that the estimates we made are 
robust. 

IV. CONCLUSIONS 

In this paper we investigated the possibility of a second-order phase transition in QCD 
in the limit of vanishing masses. When the U{1)a symmetry is restored at Tc, the finite 
temperature chiral transition, if continuous, may be described by the Lagrangian Cjj{N}) Eq- 
(1). We pointed out that the correct extrapolation at e = 1 is obtainable considering its 
U{n) X U{m) symmetric extension at fixed m. In fact the last model has a stable FF for 
n > n^{m, d) that is not accessible from the theory with n — m. The presence of a stable 
FF for light QCD with Nf = m flavors requires n~^[m, 3) > m. After showing that low order 
calculations are not conclusive, we performed a five-loop expansion that allowed to conclude 
that no continuous transition is possible for three dimensional models with U{n) x U{n) 
symmetry for n > 2. We corroborated this result with a direct three-dimensional analysis, 
namely with the pseudo-e expansion at six-loop order. 

In Ref. [4] six-loop massive zero momentum series were analyzed directly in three dimen- 
sions allowing to exclude that a FF, without any counterpart in e expansion exists (as it was 
claimed to happen for 0{n) x 0(2) models for low values of n both for > [18] and v < 
[19]). We believe that this work, together with Ref. [4], put on a robust basis the prediction 
that the transition in U{n) x U (n) models is first-order. 

In the appendix, following Ref. [4], we point out that the anomaly may lead to a con- 
tinuous transition only for Nf = 2 and large values of \wo\, instead for Nf > 3 it does not 
softens the first-order transition of the U{Nf) x U{Nf) model, for any value of Wq. 

Finally it is worth mentioning that the U{n) x U{m) models could be relevant in the 
description of some quantum phase transitions, as it happens for their 0{n) x 0{m) coun- 
terparts for Mott insulators [20]. Being n+(m, 3) > m for all m > 2, we predict a first-order 
phase transition for all those systems with m < 3, which are interesting for the condensed 
matter point of view. 
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APPENDIX A: THE EFFECT OF THE ANOMALY 

If the U {1)a symmetry is broken at Tg by the anomaly, one has to consider the Lagrangian 
(2). Since the effect of the anomaly is always (apart from the case Nf = 4) well described 
by general arguments, this appendix is very similar to part of Ref. [4] . However, we report 
such arguments here in order to make this paper self-consistent. 

In the large- A'^c limit (where Nc is the number of colors), the effect of the anomaly tends to 
be suppressed, and the Lagrangian Eq. (1) is recovered in the limit oo. The effective 

C/(l)^- symmetry breaking at finite temperature in real QCD has been investigated on the 
lattice. The U{1)a symmetry appears not to be restored at T^, but the effective breaking 
of the axial U{1)a symmetry appears substantially reduced especially above Tc (see, e.g., 
Refs. [21]). However, the Lagrangian (1) still describe a large part of the phase diagram of 
the model with broken anomaly as we show closely following Ref. [4]. 

For Nf — 2 the symmetry breaking pattern is equivalent to 0(4) — > 0(3) [2]. If the 
transition is continuous it is in the three-dimensional 0(4) universality class [2,3,22,23], 
which has been accurately studied in the literature [6,24]. Actually a continuous transition 
is expected only for large enough value of \wo\ (see Ref. [4], in particular Appendix A, 
for the phase diagram of this model). In particular the multicritical point is U{2) x U{2) 
symmetric. The phase diagram realized in light QCD may be understood only from the 
QCD Lagrangian and not from universality arguments. Lattice simulations in two flavors 
QCD favor a continuous transition consistent with the 0(4) universality class [25]. 

For Nf = 3 the determinant is cubic in $, making the Lagrangian not bounded. So the 
transition is expected to be first-order for all Wq. Lattice simulations of QCD confirm this 
expectation [26]. 

For Nf — A the determinant is quartic in leading to the three couphngs effective 
Lagrangian 

^SU{4) = ^U{4) + Wo'^ijkl'^abcd^ia^ jb^kc^ld , (Al) 

where eijki is the completely antisymmetric tensor (ei234 = 1). Such Lagrangian is not 
generalizable to an n-by-m matrix with n or m different from 4 and so we limit to consider 
the case n — m — A. The one-loop {3 functions we obtain are 

(3u{u, V, w) = —eu + 20-u^ -|- IQuv + 3f ^ -|- , 
/?u(?7,, V, w) = —ev + 8f ^ -|- 6uv — , 

(3w{u, V, w) = —ew + 6uw — 6vw . (A2) 

Such series, as in six-loop three-dimensional case [4], have no common zeros with non- 
vanishing coordinates. Anyway, differently from three dimensions, the (3^ vanishes in a 
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region of parameters, different from w = (namely the surface e — 6{u — v) =0). Higher 
loop corrections may not change the number of FP's, since they are expansions in e. 

For Nf > 5 the added anomaly term is irrelevant since it generates polynomials of degrees 
higher than four. Therefore for Nf > 5 the Lagrangians Csu{Nf) and Cu^Nj) are equivalent 
at criticality. 
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